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Abstract 

We explore the possibility of constructing stable, warped two-brane models which solve the 
hierarchy problem, with a bulk non-canonical scalar field (tachyon matter) as the source term 
in the action. Among our examples are two models-one with a warp factor (denoted as e~'^^^"^) 
which differs from that of the standard Randall-Sundrum by the addition of a quadratic piece in the 
/(cj) and another, where the warping is super-exponential. We investigate the issue of resolution 
of hierarchy and perform a stability analysis by obtaining the effective inter-brane potentials, in 
each case. Our analysis reveals that there does exist stable values of the modulus consistent 
with hierarchy resolution in both the models. Thus, these models, in which the bulk scalar field 
generates the geometry and also ensures stability, provide viable alternatives to the standard 
Randall-Sundrum two-brane scenario. 
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I. INTRODUCTION 

It is almost a decade since Randall and Sundrum (RS) hi, following earlier work 

n n n 

13|, 14, 15[, proposed a five dimensional model with a warped extra dimension, in order to 
achieve a resolution of the hierarchy problem in high energy physics and also to suggest an 
alternative to the largely ambiguous procedure of compactification, necessary in any theory 
in a spacetime with more than four dimensions. The RS two-brane model had a positive 
and a negative tension brane and we are supposed to live on the one with negative tension. 
Such negative tension branes, in general, are not stable. Moreover, it was realized that 
the modulus (representing the separation between the two 3-branes) in the RS model is 
not stable and it is easy to show that the branes would collapse onto each other due to this 
instability. In order to rescue this otherwise elegant model, Goldberger and Wise p] came up 
with a solution using bulk fields, which rendered the model stable. The visible brane tension 
however continued to be negative. Different aspects of the Goldberger- Wise mechanism have 

been discussed in numerous articles y,y,[io,[iii, 0,0,0, [0] . The stability of Horava- 



RS models) was discussed in 



Witten spacetimes (the work of Horava and Witten ^] is the string theoretic motivation for 



171]. The Goldberger- Wise (GW) resolution however assumed 



a negligible back reaction of the bulk scalar on the metric. However, it was realised that one 
could actually solve for the warp factor in the presence of the cosmological constant and the 
scalar. The back-reacted warp factor had a linear and an exponentially decaying piece and 
by adjusting parameters one could justify the exclusion of this term and work with the RS 
one. Later, it was also realised that bulk non-canonical bulk scalars (sometimes motivated 
from string theory) could rescue RS and provide stability. Examples along these lines were 
suggested and worked out in 0, Q] ■ 

In this article, we provide new examples of models with a non-canonical bulk scalar and 
without a cosmological constant. More precisely the bulk action is the tachyon matter (or 
scalar Born-Infeld) action[l^. With this assumption for the matter in the bulk, it was shown 
in 0| that one could construct a viable braneworld model. In Section II, we discuss our 
models in detail. Section III discusses the hierarchy problem and its resolution in both the 
models. In Section IV, we investigate the stability issue and obtain the inter-brane potentials 
for each model. We also make a comparison between the two models in this section. Finally, 
in Section V, we conclude with some remarks on future directions. 
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II. THE MODELS 



The braneworld models we consider in this section are, generically, two brane models. 
However, unlike the RSI two-brane scenario where there is only a bulk cosmological constant, 
we have, in our examples, bulk matter. In particular, as mentioned in the Introduction, we 
consider situations where the bulk matter is generated by a non-canonical action, namely 
that of tachyon matter. One of these models have been partially analysed in [19]. The main 
features of the model construction process is, as given below. 

We begin with the action for the bulk scalar field given by 



ST = aT j (fi^V{T)^/l + g^'^dMTdNT 



(2.1) 



where ay is an arbitrary constant, qmn being the five dimensional metric. The scalar field 
is represented by T and V(T) corresponds to its potential. The constant ut can take either 
positive or negative values. 

The full action for our model is 

S = Sg + St + Sb (2.2) 

Sg = 2M^ j ^/^Rd^i (2.3) 

- / E^O)V^^'e (2.4) 



S 



B 



where g\i} is the induced metric on the brane and r(j) is the vacuum energy density (brane 
tension) of the jth brane. There are two branes in a RS-I set-up and a single brane in RS-II 
set-up. The Sb terms incorporate the branes and the Sg term takes care of five dimensional 
gravity. In addition, there is the matter part 5*^. Variation of the total action w.r.t. qmn 
leads to the Einstein-scalar equations and that w.r.t. T gives rise to the scalar field equation. 

The five dimensional bulk metric with four dimensional Poincare invariance is taken to 
be of the form: 

ds^ = e-^f^'^^^^di^'di'' + da^ (2.5) 

Here a denotes the extra dimensional coordinate and ^ are the coordinates on a a =constant 
hypersurface. The warp factor (/) and the scalar field (T) are considered to be functions of 
a only. The extra dimension is compactified on an jZ^ orbifold and we assume the two 
branes are located at the fixed points of the orbifold, a = 0, a = Trr. With these ansatze, the 



Einstein-scalar equation without a cosmological constant reduces to the following system of 
coupled, nonlinear ordinary differential equations : 

f = (2.6) 

3 

where a — and a prime denotes derivative with respect to the coordinate a. We choose 
q;t = 1 henceforth. In general, we do not expect to obtain a general solution of the above 
equations. However, we will discuss two special solutions of these equations below. A look at 
the equation will immediately suggest that, in the bulk (i.e. excluding the brane terms) the 
relation -p- = — 2T'^ holds. Thus, postulating T one can obtain /. Further, we can multiply 
the LHS and RHS of the two equations and obtain an expression for [V(T)]^, which depends 
on /, T and their derivatives. We follow this approach while obtaining solutions. Note 
however that a /(cr) ~ |(t| would result in a constant T in the bulk (away from the branes), 
since /" would be equal to zero everywhere in the bulk. This will imply a triviality because 
the tachyon matter action would then become equivalent to having a cosmological constant 
contribution in the bulk. 

In addition to the bulk solution, one also needs to look at the equation of motion for the 
scalar field. This, for bulk tachyon matter, turns out to be: 

(l + T'2) VI + dT ^ dT ^ " ^ ' 

Integrating this equation once would give the boundary (jump) conditions on T', given as 

yn,-,^um (2.9) 



V{T) " \ dT 
where i denotes the location of the i-th brane. 

A. Model 1: 

It is possible to construct exact solutions of the above equations and boundary conditions. 
For example, one may choose: 

j{(j)^kx\a\-k\a\\^-\\ (2.10) 



Here the warp factor at cr = has a value equal to one while at a — ao — Trr its value is the 
same as that in the RS model. Notice that we have a parameter ki which, as we shall see 
below, may be assumed to depend on the modulus r. 

The full solution of the Einstein equations, i.e. T(cr) and V{T) can be found easily. These 
turn out to be: 

T{a) = -^y^ln {{h + k)ao - 2k\a\) (2.11) 

1 

V{T) = A 

The brane tensions are: 



y "0 + kao 



(2.12) 



T(i) = 24:M^{h + k) (2.13) 
T(2) = -24M^ki-k) (2.14) 

Interestingly, here both the tensions can be positive if ^ < 1 whereas if ^ > 1 one of 
them is negative. Wc shall sec later that requirements of stability and hierarchy resolution 
imply that y > 1. The tachyon field here behaves logarithmically with a and the tachyon 
potential has a complicated, but primarily exponential relation with T (which, in a sense, is 
somewhat reminiscent of the nature of the tachyon potential in the context of String theory) . 

We now need to make use of the boundary condition on the scalar field at each brane. 
Choosing the brane tensions as: 



r- e ' 

— yb arctan — + Ci 
vb 



(2.15) 



where Xi — 2.^ ^-^('^ — '^i): ^ — <^o^, ^ ~ ^ y ^ ^« integration constants, which, 
as we will see, will play a crucial role. The expression with a + sign in 2.15 is for the brane 
at (J = 0, while the one with a — sign in 2.15 is for the a = brane. With the above 
choices, we obtain the following relations by identifying the tensions given above with the 
ones obtained earher (i.e. T(i) = 24M^ {ki + k) and r(2) = -24M^ {ki - k)). 

Ci = (A;i + k)ao + ^6 arctan + j')^" (2.16) 

wb 

C2 = {ki - k)ao + ^arctan (2.17) 

Vb 

As we shall see later, the possible value of kiao (denoted as z) is dictated by the requirement 
of hierarchy resolution. The allowed value of ka^ is found from the stability requirements. 



These values can then be used to finally evaluate Ci and C2 such that the boundary con- 
ditions on the scalar field hold. Note that setting either or both of Ci, C2 to zero leads to 
inconsistencies. 



B. Model 2: 

One can easily verify that another solution exists with the following expressions for f{a) 
and T{a) 

f{a) = ^{l-e-'\^\) (2.18) 
Tia) ^ yXelH (2.19) 

with the potential V{T) given as: 
The brane tensions are: 

r(i) = 2AM^ai (2.21) 
T(2) = -24M^aie-'='^'^ (2.22) 

Notice that the brane tensions are not equal and, one of them is exponentially smaller 
(though negative in value) than the other. Also, the warp factor is decaying and has a 
value equal to 1 on the positive tension brane. Its value on the negative tension brane 
is e~'^~^^^~'^ The tachyon field grows larger as a increases. The tachyon potential, is 
nonsingular in the domain in which the tachyon field is defined. Moreover to have appropriate 
warping from Planck to TeV scale, both ai and k must be positive. 

Let us now analyse the boundary conditions on the scalar field T. Following the method 
outlined for Model 1, we first write down our choices for the brane tensions in terms of T. 
These are given as: 

(2.23) 



_^2k 
a 



^ + \nXi--ln{l + xl)+Di 



where x — ^ and Di are integration constants (the parallel of the Cj in Model 1). The 
expression with a + sign in 2.15 is for the brane at cr = 0, while the one with a — sign in 
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2.15 is for the a = ao brane. Evaluating the tensions at ai and equating them to the ones 
obtained earher in this subsection, we arrive at the following relations. 

D, = hn{2y + 1) (2.24) 
D2 = Un {2ye~^ + l) (2.25) 

where y = ^ and x = nkr. 

The first of the above relations defines Di for given y. The second defines D2 for given x 
and y. As shown later, we get y from hierarchy resolution and x from stability requirements. 



III. THE RESOLUTION OF HIERARCHY IN THESE MODELS 

We now turn towards analysing the question of hierarchy resolution in these models. 



A. Model 1 



In this model, the mass scales on the two branes are related as 

mo 

Using the definitions, z = nkir, the resolution of hierarchy requires, e" 
(or kir ~ 13, similar to the RS model). 

Furthermore, the relation between Mpi and turns out to be: 



(3.1) 



~ e-^° or 2 ~ 40 



2k p 



Erf 



Erf 



(3.2) 



where P = ^ and Erf denotes the error function. 

If z ~ 40 and p ~ 40, we find that M|,, = 0.024M^/k. Therefore assuming, M ~ Mpi 
leads to /c ~ 0.24Mp;, unlike Mpi ~ M ~ A; in the RS model (based on the relation 
Mpi ~ ^ (1 — e~^'^'^^''=) in RS and for moderate fcr^ ~ 12). Note however, that the k in 
our model is not the same as the RS 'k'-which is attached to the linear part in the RS 
warp factor. In our model, the parallel of the RS 'k' is ki. Since ki = kp ^ 40k we get 
ki ~ 0.96Mpi. Thus, one can set fci, M and Mpi of the order of Mpi while k needs to be 
set at an order of 10~^Mp;. Therefore, in this model, an extra small hierarchy, somewhat 
similar to the ADD scenario exists lal and cannot be avoided. 



B. Model 2 



We note that the relation between mass scales on the branes at a = and o" = vrr is 

m = m^e-f^^'^ = e-s^W(i-^"^) (3.3) 
To solve the hierarchy problem one needs — ~ e~^°. Thus, ^ must be around 40. 

mo ' k 

In addition, the relation between Mp^ and turns out to be: 

Ml, = -^e-'y [Ez[2y] - Ei{2ye-% (3.4) 

where Ei denotes the exponential integral function. 

With y = 40 and x moderate (about 3 or larger, as we will see from the stability analysis), 
we obtain, for Mpi ~ M, k ~ 0.025Mp;. As in the Model 1, we have the parameter a\ as 
the parallel of the RS 'k' and this takes on the value a\ ~ 40A; ~ Mpi. Similar to Model 1, 
an extra small hierarchy seems to exist here too (through the value of the k in our model) 
and the resolution of hierarchy seems to have an ADD-like flavour 



IV. STABILITY ANALYSIS: INTER-BRANE POTENTIALS 

We shall now investigate the stability of our two-brane models. 

Let us first look at the original RS-I model (with only a bulk negative cosmological 
constant) and construct the interbrane potential. This turns out to be (omitting the / d^^) 

V{x) = -ISM^k (1 - e"^^) (4.1) 

where x = nkr and the warp factor is just e"^'*^'. Note that there is only one free parameter 
here (i.e. x) and the potential does not have any minimum. 

This potential will be modified, if we also take into account the five-dimensional gravita- 
tional part of the action and integrate it over the extra dimension. In summary, the gravity 
part is given as J y/g^Rd^C, which has a piece which would give 4D gravity upon integrating 
over extra dimensions (provided we have a curved brane) and another piece (proportional 
to / e-^f (8/" - 20/'^) dad^O- ^ is a debatable issue as to whether we should include this 
piece in the potential or not. We, however remain non-commital on this and discuss both 
scenarios (i.e. with and without the gravitational piece). Including this so-called gravity 
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part into the inter-brane potential for RS-I, we find that 

V{x) = AM^k (1 - e~^^) (4.2) 

The potential is now positive definite but, as before, has no minimum. 

The above analysis shows that, at least for RS-I, by including the gravity part one cannot 
make the branes stable. Goldberger and Wise Gj therefore included an extra bulk scalar 
field which was crucial in order to achieve stability. 

In our models, we have bulk tachyon matter which generates the spacetime geometry. In 
both the models discussed above we have two parameters( z, p in Model 1 and x, y in Model 
2). Following GW, we consider the matter action (tachyon-matter) + the two brane actions 
and insert the solutions for T{a) and /(cr) in them. Thereafter, we integrate over the extra 
dimension to obtain the effective inter-brane potential as a function of r. 

As mentioned before, in either model, we have two independent parameters and, further, 
a pair of constants which depend on these parameters. We adopt the following strategy to 
find the allowed values of the parameters and the constants. Let us assume that in either 
model the parameters are related functionally. We do not provide any compelling reason 
for doing so-however there is no reason either why such a relationship may be disallowed. 
In other words, in Model 1, we have p{z) and in Model 2, y{x). We can now choose simple 
forms of these functions (preferably trigonometric) in order to have a minimum in either 
potential. Let us now focus on the two models separately. 



A. Model 1: 



Firstly, let us follow the procedure outlined above to obtain the potentials. Define p = 
ki/k as before, so that 

/(a) = k 

Now let z = nkir and p = p{z). Then, the expressions for the potentials, including and ex- 
cluding the Ricci scalar contribution in the action are, respectively (omitting the integration 
over d'^x), 



p\a\ 



(4.3) 
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Without gravity term contribution: 



V{z) = -24M' 




+ 




(4.4) 



Jo 



With gravity term contribution: 



Vrnodiz) = 8M^kip+l)-8M^k{p-l)e-^'- [ 



Jo 




1 1 

p 




(4.5) 



where, in the above integrals x — ^. Note, both the integrals above can be done-we do 
not write them out here because the expressions are complicated and it is more useful to 
plot the potentials directly. It should also be realised that the dominant contributions to 
the potentials essentially come from the first two terms in each of them. 

The strategy to find possible stable values of z will be as follows. Note that the hierarchy 
resolution requires z — 40. Let us choose p{z) = pocos^ z. This restricts p to have values 
between and po- Further it ensures the presence of minima in V{z) since the linear in p 
terms will contribute to their existence. Using this form of p{z) (with po — 40, say) in either 
of the above expressions for the potentials we find that there is a minimum around z — 40 
(the value required for hierarchy resolution). We show the potentials in Fig. 1 and Fig. 2. 
Since z and p values are now known and we can easily find out Ci and C2 by substituting 
these values in the expressions for the Cj given earlier. 

B. Model 2: 

For the analysis of the stability of branes in Model 2, recall, the dimensionless variables 
y = ^ and X = irkr. Omitting, as before, the explicit writing of the integration over d'^^ we 
obtain the following effective potentials: 
Without gravity term contribution: 




(4.6) 



With gravity term contribution: 
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V(z)/M"3 k 




FIG. 1: Potential V{z)/M^k vs. z (without gravity term contribution) for Model 1. Note a 
minimum at z ^ 41. 

V_mod (z)/M"3 k 




32 34 36 38 40 42 44 



FIG. 2: Potential Vmod{z)/M^k vs. z (with gravity term contribution) for Model 1. Note a 
minimum around z ~ 39.5. 



1 - e-^^e-^e^^"" 



(4.7) 



The dominant contribution in each potential comes from the linear in y term. In either 
expression, we now choose y{x) — ygcos^ x (with yo ~ 40). Once again, this choice restricts 
y to values between and 40. With this choice, we now plot the potentials and see whether 
there are minima at finite values of x. 
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The figures (Fig. 3 and Fig. 4) below demonstrate that we do indeed have minima in 
the potentials such that y is close to 40 at some finite positive x. 

V (x) /M^3 k 




FIG. 3: Plot for V{x) /M^k vs. x (without gravity term contribution) for Model 2. Note a minimum 
around x ~ 3. 

V_mod (x)/M"3 k 




FIG. 4: Plot of Vmod{x)/M^k vs. x (with gravity term contribution) for Model 2. Note a minimum 
around x ~ 4.7. 

Thus, as in Model 1, knowing x and y we can easily find out Di and D2 from the 

abovementioned expressions. 

C. Comparison between the two models 

It is reasonable to ask-which among the above two is a better model. A very clear answer 
to this question does not seem to exist, at this stage. We list below a few points which might 
shed hght on this aspect. 
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• The second model is very different from RS in terms of the functional form of the warp 
factor. The first one can, in a way, be thought of as a perturbation. 

• At the level of solutions (without invoking any constraints that arise from hierarchy 
resolution or stability). Model 1 seems to allow the existence of two positive tension branes, 
while Model 2 does not. This feature of Model 1 differs from what we obtain in standard RS 
where there is no scope of having two positive tension branes even at the level of solutions. 

• In the first model the scaling of masses occurs in a way different (super-exponential) 
from RS whereas in the second model the scaling is the same (exponential) as for RS. 

• The nature of the tachyon potential in the first model is far removed from whatever 
string theorists write down from varied considerations (which, in some sense may not be 
applicable here because there arguments are based on cosmological considerations-rolling 
tachyon etc.). In the second model, by virtue of being exponential (not exactly, but in some 
way related) the tachyon potential is perhaps a bit closer to the stringy ones. 

• The stability issue is solved in both the models by incorporating a similar logic (i.e. a 
functional realtionship between the two dimensionless parameters). Neither model can be 
called as preferred, from this perspective. 

• For hierarchy resolution, we have seen that Model 2 is different in nature than Model 
1, which is closer to RS. This is because in Model 1, the value of p does not affect hierarchy 
resolution as long as the choice of p{z) results in a minimum in V{z) at the right value of 
z (required for the solution of the hierarchy problem). On the other hand, in Model 2, a 
proper value of y is required for hierarchy resolution and the function y{x) should attain 
that value at a minimum in the potential. 

• Further, in the resolution of hierarchy in either model we seem to be introducing an 
element of hierarchy through the values of the extra parameters (the k in Models 1 and 2). 
This aspect seems to be closer to the hierarchy resolution in the ADD model jsj , though it is 
not, in any way dependent on the number of extra dimensions, as in ADD where it happens 
through the factor associated with the volume of the extra dimensions. 

In order to facilitate a comparison we tabulate the essential results of Models 1 and 2 in 
Table I. 
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Feature 


Model 1 


Model 2 


Warp factor 


f{a)^h\a\-k\a\ - l) 


/(a) = f (1 - e-^l'^l) 


Taeliyoii fic4d 


^(^) = -| yf ln((A-i + /,-)ao - 2k\a\) 




Tachyon Potential 


i 

y(T) = A e V-o' +A;(To e V-o' 




PaTSTn ptprs 


k fci r 


0,1 k r 
'-''1) ^) ^ 


Rranp tpnsions 


Mil ^TiJ. vj. V 'I'l 1 ' / 7 

r/,^ = -24M^ ffci - k) 


Tn', = 24M^ai 
rr9^ = -24M^aie"'""' 


Reduced parameters 


z = Tr/cir, p = ^ 


y = X = nkr 


B.C. and constants 


Ci,2 (Eq. 2.16, 2.17) 


L>i,2 (Eq. 2.24, 2.25) 




9 

T) — T)r> pn^ 7 


2 


Hierarchy resolution 


?D rsj TIpV T77n A/fru 7 r\J 40 


T ~ 3 to 5 


Mpi, M relation 


Eq. 3.2 gives for 
M ~ M d; ~ A;i A; ~ 024Md; 


Eq. 3.4 gives for 
M ~ Md; ~ qi A; ~ 025Md; 


Inter-brane potentials 


Eq. 4.4, Eq. 4.5 


Eq. 4.6, Eq. 4.7 


Stability 


Pot. minimum around 2; ~ 40 

p = Po ~ 40, cos^ 2; ~ 1 


Pot. minimum around x ~ 3 to 5 

y = Ho ^ 40, cos^ x ~ 1 


Deficiency 


Additional small hierarchy in k 


Additional small hierarchy in k 



Table 1: Summary and comparison of the two models. 



V. REMARKS AND CONCLUSION 

In this article, we have tried to build two-brane models with bulk matter, which solve 
the stability and hierarchy problems. 

How different are our models from the standard Randall-Sundrum ? From Table 1, we 
can write down the warp factors for the two models using the values of the parameters 
quoted there. These turn out to be: 
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Model 1: 



f{a) = Mpi |a|-0.024|a| i-^Mp,-l 




) 



(5.1) 



Model 2: 



fia) = 40 (1 



^-0.025Mpi\a\ 



(5.2) 



If we plot the above two warp factors along with the linear-in-a RS one we will see that the 
deviations are very small and the overall decaying nature is maintained. So, where does this 
small difference make a crucial contribution? To see this let, us go back to the expressions 
for the inter-brane potentials. In either model, the inter-brane potentials are dominated by 
the hnear terms {z in Model 1 and y in Model 2), both of which are essentially the ratio of 
each of two parameters {ki, k in Model 1 and ai, k in Model 2). The existence of these linear 
terms therefore depends on the existence of more than one parameter (apart from r) in either 
model, which, in turn, is made possible through the presence of bulk matter. Postulating 
a relationship between the parameters (i.e. y{x) or p{z)) thus becomes a possibility (unlike 
the RS case) and hence shows a way of achieving stabihty. 

The tachyon matter scalar field is responsible for providing the source for the bulk metric 
and keeping the branes stable. The crucial point in our analysis has been the fact that among 
the two free parameters in each of our models (x,y and z,p), we have invoked relationships 
(i.e. y{x), p{z)). This has facilitated the resolution of the stability and hierarchy problems 
simultaneously. We agree that we do not have a good reason behind choosing y{x) or p{z) 
though, as shown, trigonometric functions seems to be the preferred choice in order to have 
minima in the potentials and hence, stability. 

The issue of whether to include the contribution of the gravitational part of the action 
in defining the effective inter-brane potential remains an open issue. We have preferred to 
be non-commital on this in our paper and that is why we include results for both scenarios. 
However, we may note that we do not see any major difference in our results after including 
the gravitational part of the action. The only difference seems to be the fact the in the 
presence of the gravity part the potential in all cases is positive, while, otherwise, it is 
negative. If we had seen a more significant difference (i.e. say, the existence of minima 
with the gravity part and no minima without or vice versa), we might have had a case for 
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choosing one over the other. At this stage, we do not make any definite statement on this 
aspect. 

We conclude this article by stating that our analysis of stability and our proposal on a 
functional relationship between the parameters in our models, surely needs more justification 
at a more basic level. In addition, it is crucial to address field localisation, massive modes, 
KK graviton effects, corrections to Newtonian gravity and other phenomenological issues 
using the new warp factors introduced in this article. We hope to address some of these 
issues in our future investigations. 
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